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1. Coulomb gauge at leading order 

The first part of this talk concerns the construction of a leading order truncation to the Dyson- 
Schwinger equations of Coulomb gauge QCD [|T]]. Let us begin by considering the following 
(standard) functional integral 

' 2^e''^\ y = 9^ + J dx{E^-B^)/2, (1.1) 

where the action (^5^) is split into a quark component, ^5^^, and the Yang-Mills part. i^<I> generically 
denotes the integration measure. The chromomagnetic field, B, will not concern us in the following. 
The chromoelectric field, E, is given by (superscript indices a,... denote the color index in the 
adjoint representation) 

E" = -doA"-D"''A^^, b"^ = 5"''W - gf'^A' , (1.2) 

where D is the spatial component of the covariant derivative in the adjoint color representation (the 
/ are the usual SU {Nc) structure constants). We work in Coulomb gauge (V • A = 0), for which 
the corresponding Faddeev-Popov (FP) operator is — V • D. There are two important points: the 
FP operator involves purely spatial operators and the chromoelectric field is linear in the temporal 
component of the gauge field, Aq. We now convert to the first order formalism by introducing an 
auxiliary field n via the identity 

exp{,/..v.V2}^/«exp{,/.4-.V2-*«.^i}. 

The field n is split into transverse iW -n^ = 0, henceforth we drop the _L) and longitudinal (V0) 
parts. Since the action is now linear in Aq, we can integrate it out, to give 

Z=y"^<I>5(^V-A)5(^V-7r)Det -V-D 5 (^V -D^ +p) e''^', p" = gf^^'A'' -tz' + gq[fT'']q, 

(1.4) 

where p is the color charge (including the quark component, with quark field q and the Hermitian 
color generator T"). The ^-field can also be integrated out to cancel the FP determinant (Coulomb 
gauge is formally ghost free). However, noting the temporal zero modes of the FP operator, i.e., 
those spatially independent fields for which —V • {xq) = 0, one is left with [§] 



Z = j^d (v • a) 5 (v • tt) 5 (^y^ dx^ e'-^\ y" ^ jdx 



papahphi2 



(1.5) 



In the above, one sees that there are two transverse degrees of freedom for the gluon and the total 
color charge is conserved and vanishing. The Coulomb kernel F = [— V-D]^'(— V^)[— V-D]^' is 
nonlocal in A, so we make the leading order truncation whereby it is replaced by its expectation 
value, which is related to the temporal component of the gluon propagator: F -^<F>^ Wbo. It is 
known that in Coulomb gauge on the lattice. Woo is infrared (IR) enhanced, going like o /q^ but 
with a coefficient a somewhat larger than the Wilson string tension (see e.g., Refs. [^, Q]). The 
charge conservation term is rewritten in the limiting form of a Gaussian, mimicking the Coulomb 
term: 

5(^1 dxp^ ~ lim ^exp ^-1/2 J dxdyp;^d'''''^6{xQ-yo)py^, (1.6) 
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and such that we now have instantaneous four-point interaction terms including Taajui and F; 



qqqq- 



dxdy 



p"J"'F„p^/2 , g'CFF{q) = {2Ky^d{q) + %na/q 



7?* 



(1.7) 



{Cp = (A'^^ — \ )/2Nc)- This interaction contains the charge constraint and leads directly to a linear 
rising potential with a string tension a. To complete the leading order truncation scheme, we 
restrict to one-loop terms in the following equations and disregard all but the F interaction terms. 



2. Gluon gap equation 

In the first order formalism, the transverse spatial gluon degrees of freedom (A, k) have been 
separated such that there are three propagators Waai], ^mtij and WAnij are the spatial indices), 
correspondingly with three proper functions r^,^, Tjijiij and TAnij related via a matrix inversion 
structure (see, e.g., Ref. [||]). Since the interaction content of our truncated system is instantaneous, 
the energy dependence of these functions is trivial (and the mixed functions will play no role in 
the discussion here). There are two scalar dressing functions of interest, both functions of spatial 
momentum: FaaC^ ) and rj[j[(P-). The spatial gluon propagator Waa has the form {Wji„ is similar) 



WAAijik) = itij(k) 



[kl-k2rAA{k^)r^^{k^) + iO, 



(2.1) 



{tij is the transverse spatial momentum projector). The truncated Dyson-Schwinger equations have 
the mnemonic form (omitting kinematical factors etc.) 

r„nif) ~ 1 + / dkF{p-k)WAA{k), Taa{p^) ~ 1 + / dkF{p-k)W^^{k). (2.2) 

The charge constraint term of the interaction F (i.e., the term ^ 'io8{p — k)) immediately tells us 
that both and T^j^ are divergent as ^ — ;> oo (there is also an IR divergence), meaning that the 
gluon self-energy is infinite and the propagator poles are shifted to infinity. This has the natu- 
ral interpretation that one requires infinite energy to create a (colored) gluon from the (colorless) 
vacuum. If, however, we consider the static gluon propagator , written as 



^mtik^^ 



^WAAiiik)=tij{k)^, Gk 
2k ^''^ 2\k\ FAAik^) 



(2.3) 



then we can combine the Dyson-Schwinger equations to get the gluon gap equation 



Gn 



1 + ' 



g^Nc fdkF{p-k) 
{2ny\k\ 



tji{p)tij{k) 



k^ G, 



(2.4) 



This equation has previously been derived in the Coulomb gauge Hamiltonian approach [^. The 
dressing function for the static propagator, G, is IR finite and independent of the charge constraint. 
Solving numerically (in units of a), one sees that the solution has the form Gx = x/(x+ K^) for 
an IR constant 'mass' function k{x) and where x = Ic. is logarithmically dependent on the 
numerical ultraviolet (UV) cutoff A (dimensions of [mass]^), despite the fact that the interaction 
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Figure 1: [left panel] Kx as a function of x — l? and [right panel] ^{x') as a function of x' —x/a for various 
values of the UV-cutoff A. All dimensionfull quantities are in appropriate units of the string tension, a. See 
text for details. 



has the form l/q^. K is plotted in the left panel of Fig. ^ However, defining a = k{x = 0) and 
introducing the scaled variable x' = x/a, one finds that k{x!) = k{x = x'a) — a is independent of 
A, shown in the right panel of Fig. |T} It turns out in general that by simply writing all dimension- 
full quantities in units of the (dynamically generated) gluon mass function at some point, one may 
construct A-independent dressing functions (and subsequently G) without introducing a renormal- 
ization constant 



3. Quark gap equation 

Given that the interaction content arising from the Coulomb term couples to the gluon and 
quark charges in the same manner, the quark sector turns out to be very similar to the gluon sector 
within the leading order truncation scheme considered here. The instantaneous character of the in- 
teraction leads immediately to the following form for the quark propagator in terms of two dressing 
functions A and B (both functions ofk ): 

[kl-k^Al-Bl + iO+] 



A possible term ~ 'fk^y-k does not appear, just as in the perturbative case [^. The static quark 



propagator, W^'"' , can be written in terms of the mass function, M, and quasiparticle energy, ft): 



= I '^W,,(k) = f±^, M, = ^, col = +Ml (3.2) 
The Dyson-Schwinger equations for the dressing functions A and B have the mnemonic form 

Apr^ 1 + J dkF{p-k)/(Ok, Bpr^m + J dkF{p-k)Mk/(Ok, (3.3) 
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Figure 2: [left panel] Quark mass function, M{x), and [right panel] dressing, M{x) — m, plotted as functions 
of X — p for a range of quark masses. All dimensionfull quantities are in appropriate units of the string 
tension, CT. See text for details. 



showing us via the charge constraint that the quark self-energy is divergent (like for the gluon) and 
one requires infinite energy to extract a single quark from the vacuum. However, combining the 
equations in terms of the mass function, M, leads to the Adler-Davis gap equation ^ 



1 



m- 



dkF{p — k) 



(34) 



The mass function is IR finite and independent of the charge constraint. While the above equation 
was originally derived for chiral quarks (in the Hamiltonian formalism), in the leading order trun- 
cation scheme presented here one can show [jl]] that it also reproduces the Coulomb gauge heavy 
quark limit (in the absence of pure Yang-Mills corrections) [^]. The mass function is plotted on the 
left panel of Fig. ^. One can see that chiral symmetry is dynamically broken, although the chiral 
condensate is too small [|l^] (this can be improved by considering the spatial quark-gluon vertex 
[11]). In the right panel of Fig. ^ the dressing M(x) — m is plotted. As the quark mass increases 
the dressing initially also increases, but for heavier quarks becomes smaller and in the heavy quark 
Umit, M 



m. 



4. Bethe-Salpeter equation 

Within this leading order truncation scheme, it is possible to study the quark-antiquark Bethe- 
Salpeter equation for (color singlet, flavor nonsinglet) pseudoscalar and vector mesons with arbi- 



trary quark masses [|10|j. The pseudoscalar case will be discussed here - the vector case is similar. 
In the Coulomb gauge rest frame, the Bethe-Salpeter vertex for pseudoscalar meson can be written 
(omitting flavor factors) 



5 



Gap and Bethe-Salpeter equations in Coulomb gauge 



Peter Watson 



where Pq = Mj,^ is the total energy squared (at resonance) for the quark-antiquark pair and p the 
spatial momentum flowing along the quark line. The dressing functions F, all have the argument 
p . There are two basic quantities of initial interest (trace over Dirac matrices): 




Trj 



dk \ 'fP{)f 



{2n) 



W+{k+)rps{k;Po)W^-{k 



qq\ 



(4.2) 



where k'^ represents the energy and spatial momentum argument {k^ ±^0/2,^) and the two quark 
propagators W-^^ correspond to bare quark masses ni^. fps is the pseudoscalar meson leptonic decay 
constant, hps is related to fps via the axialvector Ward-Takahashi identity (AXWTI) [10, and 
this can be compared to the Gell-Mann-Oakes-Renner relation in the chiral Umit 



Mlsfps={m^ + m )hps, hps' 



<qq> If PS, 



(4.3) 



indicating that hps is a generalization of the chiral condensate to finite, arbitrary mass quarks. 



Evaluating the trace and energy integrals for the right-hand side of Eq. {\.2), one obtains spatial 
integrals over a combination of terms involving IR divergent quantities such as A. However, fps 
and hps must be IR finite. Assuming the form 



fps = 2iNc 



dk 



PS ■ 



■llNr 



dk 



-[(0+ + (0,^]hk, (4.4) 



the combinations of divergent factors are then contained within the two functions / and h. Here 
is where Coulomb gauge does something special: when one expands the truncated Bethe-Salpeter 
equation 



Tps{p;P'' 



dk 
(2^ 



F{p-k)fW^+{k+)rps{t,Po)W^-{k-)f, 



(4.5) 



the right-hand side takes the mnemonic form F,- ~ J F[...] [//. or hj^] where the terms represented by 
the dots ([. . .]) involve combinations of only the finite functions or (O^. The Bethe-Salpeter 
equation can thus be rewritten in terms of only / and h. The equal mass case is 



hr 



1 



p 



dkF{p-k) 



h^ — h 



p-k 



fp = hp + -g Cf 



1 2 fdkF{p-k) \ [p-k + MpMk 



{27iy (Ok 



fk 




(4.6) 



The arbitrary mass case has a similar form. The corresponding vector meson equation is also sim- 
ilar, but involves four functions. One can see that the above form for the Bethe-Salpeter equation 
thus behaves like the previously discussed gap equations for G and M, where the charge constraint 
and IR divergences cancel and despite the interaction, the functions / and h are finite. The equa- 
tions can be compared to those of, for example, Refs.[p3|]. 

Turning to the results, the normalized (see Ref. [|0[) pseudoscalar and vector meson dressing 
functions are plotted for the chiral quark case in Fig. || and one sees that indeed, the functions 
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Figure 3: [left panel] Pseudoscalar and [right panel] vector normalized vertex functions with (equal) chiral 
quarks, plotted as a function of x = P. All dimensionfuU quantities are in appropriate units of the string 
tension, a. See text for details. 
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Figure 4: Pseudoscalar and vector meson masses [left panel] and leptonic decay constants [right panel] with 
equal mass quarks, plotted as a function of the quark mass. All dimensionfuU quantities are in appropriate 
units of the string tension, a. See text for details. 



are all finite. In Fig.^ the meson masses and leptonic decay constants for equal quark masses are 
plotted (in units of the string tension a). Inserting typical values for a [|l0||, it becomes obvious 
that whilst dynamical chiral symmetry breaking is visible {Mps ~ -Jni as m — )• 0), the leptonic 
decay constants are too small, as is the mass-splitting between states for larger quark masses. 
Looking at the case for one fixed chiral quark plotted in Fig. ^ one sees that both the pattern 
for chiral symmetry breaking {Mps ~ ^/fn as m — )• 0) and the leading order heavy quark limit 
if PS V ^PS ) fv V My ~ const, as m — )• oo) are present. The leading order Coulomb gauge truncation 
scheme thus qualitatively accommodates both chiral and heavy quark physics. 
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Figure 5: Pseudoscalar and vector meson masses [left panel] and fps\/Mps, fwJMy [right panel] with 
one fixed chiral quark, plotted as a function of the other quark mass. All dimensionfull quantities are in 
appropriate units of the string tension, a. See text for details. 
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